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The local density of states of Ag(lll) has been probed in detail on disordered terraces of varying 
width by dl/dV-mapping with a scanning tunneling microscope at low temperatures. Apparent 
shifts of the bottom of the surface-state band edge from terrace induced confinement are observed. 
Disordered terraces show interesting contrast reversals in the dl/dV maps as a function of tip-sample 
voltage polarity with details that depend on the average width of the terrace and the particular edge 
profile. In contrast to perfect terraces with straight edges, standing wave patterns are observed 
parallel to the step edges, i.e. in the non-confined direction. Scattering calculations based on the 
Ag(lll) surface states reproduce these spatial oscillations and all the qualitative features of the 
standing wave patterns, including the polarity-dependent contrast reversals. 
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I. INTRODUCTION 

The fcc(lll) surfaces of noble metals exhibit a large 
sp-band gap in the projected bulk band structure along 
the T-L linei^ These gaps reach down below the Fermi 
energy at the center of the surface Brillouin zone and 
support an occupied free electron surface-statei&i Elec- 
trons occupying these Shockley-type surface-states form 
a two-dimensional nearly free electron gas parallel to the 
surfacei^iSiL Electrons are confined to the vicinity of the 
top layer by the vacuum barrier on one side and a band 
gap in the bulk states on the other side. 

Scanning tunneling microscopy (STM) and spec- 
troscopy (STS) are particularly sensitive tools to study 
these surface-states that dominate in the zone be- 
tween the tip and the surfaced These techniques 
map the local electronic density of states (LDOS) in 
real space. STM has provided real space observa- 
tion of electrons in surface-states and of their interac- 
tions with adsorbateSfSiiS steps j2ii2iiiii2ii&i4, and other 
structures2ii2ii5i±& in the form of standing wave patterns. 
These spatial oscillations are quantum mechanical inter- 
ference patterns caused by the scattering of electrons 
in the two-dimensional electron gas from defects. The 
standing wave patterns contain information about the 
surface-state dispersion and they give insight into the in- 
teraction between surface-state electrons and scattering 
sites on the surface AiSi 

STS images of the local differential tip-surface con- 
ductance (called dl/dV maps) have been used to in- 
vestigate how the surface-state interacts with a single 
straight atomic stepiii^ and with a pair of adjacent 
straight stepsiii*i£ From the measurement of spatial os- 
cillation periods normal to the step it was possible to 
extract Fermi wave vectors^ and complete energy disper- 
sion curvesi^S Thus, the surface-state on clean, disorder- 
free systems is well characterized. However, growth and 
catalysis, which are both influenced by the electronic 



structure of the surfacej2ii22i2ii also occur on surfaces 
where disorder is present. In particular, the LDOS mod- 
ulations on irregular step arrays^* can have a dramatic 
impact on these processes. Therefore, it is important 
to understand and to be able to predict the electronic 
structure on disordered terraces of varying width limited 
by steps with irregular edges. A study of the electronic 
structure of such terraces is the focus of this paper. 

In a recent publication, we have investigated the de- 
pendence of the LDOS near the Fermi-energy, ~Ep, on the 
width of terraces on Ag(lll) at 7 KS We showed that 
the apparent surface-state band minimum shifts mono- 
tonically towards the bulk Ep with decreasing terrace 
width leading to a depopulation at a mean terrace width 
L = 3.2 nm, in quantitative agreement with a Fabry- 
Perot model with energy dependent finite asymmetric 
reflectivity of the step edges. In addition, a switch 
from confinement by terrace modulation to step mod- 
ulation was observed at \f/2 (L=3.7 nm), half of the 
Fermi wavelength on an infinite terrace. This switch 
occurs on Cu(lll) surfaces also at \p/2 as proven by 
photoemission 

In this paper, we present STS data revealing spatial 
oscillations of the local density of states on Ag(lll) for 
narrow disordered terraces of varying width up to 8.5 
nm. Apparent band ed ge s hifts of the surface-state on- 
set (as reported in Ref. |25j) lead to conductance varia- 
tions on terraces of different width at the same energy. A 
standing wave pattern also evolves parallel to the steps, 
i.e. in the non-confined direction. Its wavelength follows 
the dispersion relation of an electron in the clean, infi- 
nite surface-state. In order to understand the wave pat- 
tern, theoretical simulations based on scattering theory 
are performed for disordered terraces. These calculations 
successfully reproduce the observed oscillations parallel 
to the terraces as well as interesting tip-sample voltage 
polarity dependent effects. 

This paper is organized in the following way. In Sec.lITI 
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we discuss the details of the experiment, including sam- 
ple preparation, the types of measurements made, and 
the main features observed in the data. In Sec. IIIII wc 
discuss the theoretical model used to interpret the ex- 
periments and its assumptions. We also present the re- 
sults of model calculations and discuss their connection 
with experiment. Finally, in Sec. IIVI we give the main 
conclusions of our work. 



II. EXPERIMENT 
A. Sample preparation and data acquisition 

The experiments have been performed in ultrahigh vac- 
uum with a low temperature scanning tunneling micro- 
scope that operates at temperatures between 6 and 300 
K. The single crystalline Ag(lll) surface is cleaned by 
sputtering and annealing cycles. The surface separates 
into large flat terraces and step bunches with an average 
terrace width of 4 nmiSi In this paper, we investigate the 
spatial oscillations of the LDOS in these stepped regions. 
Measurements are performed at 7 K. 

The topographic images are taken in constant current 
mode. dl/dV spectra are recorded with the lock-in tech- 
nique; the ac tunnel current is driven by a 4 mV signal 
added to the junction bias with a frequency of 327.9 Hz, 
381.8 Hz, or 738.1 Hz. Voltages are applied to the sample 
with respect to the tip. Thus, for a negative voltage the 
occupied side of the spectrum is probed. dl/dV maps are 
recorded simultaneously with the topographic images at 
the same voltage by recording the lock-in signal at each 
pixel. It is not possible to run the STM in constant height 
mode in imaging a step array of more than 1 nm change in 
height. Thus, the constant-current image recorded simul- 
taneously with the dl/dV map is used for an adjustment 
in tip-sample separation. This leads to a superposition 
of the dl/dV signal with a variance in height at the step 
edges. 

A dl/dV map of 256x256 pixels is recorded in about 80 
minutes. Series of dl/dV maps therefore requires a high 
stability STM, which our system provides. Typically, 
the drift is less than 0.3 nm during a series of 10 dl/dV 
maps, i.e. 13 hours. Major drift problems arise from the 
increasing temperature at the end of helium evaporation. 
Therefore, measurement series on the same spot of the 
crystal are limited by the helium amount available in the 
bath cryostat. 



B. Experimental results 

For a defect-free terrace of ideally parallel steps, the 
surface-state electrons on each terrace are confined per- 
pendicular to steps but are free parallel to them. Thus, 
the wave functions of surface-state electrons are separa- 
ble, with standing waves normal to the steps and plane 



waves parallel to the steps. This leads to a wave pat- 
tern perpendicular to the steps. Fig. ^ shows these well 
known spatial oscillation of the surface-state electrons 
perpendicular to pre-existing step edges on a large terrace 
of 51 nm in width. In addition, surface impurities act as 
point defects and lead to radial wave patterns centered 
on the impurities. At 32 mV the waves have the expected 
wavelength of A/2 = 3.5 nm (See Ref. [U). At 7 K five 
maxima are discernable already in the topographic im- 
age. Note the irregularities of the plane waves introduced 
by the natural deviation of the step edge position from a 
straight line. 

The LDOS (as recorded in dl/dV spectra) displays a 
step-like function marking the sharp increase of differ- 
ential conductivity around the surface-state band onset 
(Fig.HJj). Following Ref. 0, we measure the onset of 
the surface-state, U on = {Utop + Ubot)/2, and a width, 
A = Utop — Ubot, of the onset by continuing the slope at 
the middle of the rise to the bottom and the top of the 
rise as indicated in Fig.QJ). The values of U on — — (68±4) 
meV and A = 8 meV are in good agreement with earlier 
measurements i& 

Confinement on narrow terraces leads to both a shift 
in energy of the surface-state related peak and a broad- 
ening of itiS^ Moreover, additional maxima develop. The 
example shown in Fig. ^i, taken on a terrace of 6.1 nm 
width, displays onsets at -53 meV and 32 meV with onset 
widths of A_53 me y = 23 meV and A^meV = 32 meV, 
respectively. These two peaks correspond to the first two 
transverse modes of the surface state confined by the ter- 
race edges. 

The spectra vary for terraces of different width, but 
also for different positions on the same terrace (Fig. |3J) 
with slight variations in apparent surface state onset en- 
ergy and large variations in dl/dV intensity (Figs. 03 and 
c). In particular, close to the step edges, Ubot shifts to 
lower values influencing the onset value U on (Fig. |5Jd). 
We have shown before that these variations close to the 
step edges are in accordance with a Fabry-Perot model 
with the steps acting as finite, energy-dependent barriers, 
so this is no indication of a shift in energy^ As a conse- 
quence of strong scattering of the surface-state electrons 
by the step potential, the amplitude of the surface-state 
is smaller close to the step edges. 28 Thus, the peak inten- 
sity of dl/dV varies for the same terrace width depending 
on the distance from the step edges (Fig.[2t), being high- 
est in the middle between the step edges for the first peak 
in dl/dV and showing two maxima for the second peak 
in dl/dV on the largest terrace. Only the widest terrace, 
the second one in Fig. shows two peaks in dl/dV. 

The second peak arising in dl/dV due to local con- 
finement on a single terrace should show a quadratic de- 
pendence of energy with terrace width according to a 
particle-in-a-box model for the simple case of infinitely 
high potential wallst^i E gnc i — h 2 Tt 2 /2m e ffL 2 , with m e ff 
the effective mass and L the terrace (confinement) width. 
The plot of the energy of the second peak relative to 
apparent surface-state onset vs. 1/L 2 (Fig. [3J shows 
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indeed a linear dependence for terrace widths > Af/2, 
i.e. for those terraces where the Fabry-Perot model is 
applicable. 25 From the prefactor of (6.371 ± 1.794) eV 
nm 2 , we determine m^TJ = (0.59 ± 0.17) m e , which de- 
viates from the effective mass of electrons in the surface- 
state band where m e // = 0Am e . 

To investigate spatial variations of the differential con- 
ductance in more detail, we have recorded dl/dV maps 
on these step arrays at different energies. Some examples 
are displayed in Figs. 0]to[7| 

Below the onset of the surface-state, dl/dV maps dis- 
play bright lines and spots on a continuous background 
(Fig. to d) . The bright lines correspond to the surface 
steps and are an artifact of the constant current scanning 
process (see above). The brightest spot in the dl/dV 
maps corresponds to an impurity at the step edge. It is 
visible in the topography as a protrusion imaged 20 pm 
higher than the surrounding step edge (Fig.0^,). This im- 
purity can be used as a reference point to compensate for 
thermal drift. The two less bright spots are not visible in 
the topographic STM image in the contrast shown. Con- 
trast enhancement shows them as protrusions of only 4 
pm in height. We therefore attribute them to subsurface 
impurities. 

With increasing voltage a difference in contrast devel- 
ops between terraces of different width. For the broadest 
terrace "3" , there is a contrast switch from dark to bright 
between -50 mV and -40 mV (Fig. 0^ to h). Narrower 
terraces switch contrast at voltages closer to Ep. The 
dependence of the switch in intensity on terrace width 
results from the different apparent surface-state onsets 
(see Fig.|2Ji). On the largest terrace "3" a wave pattern 
evolves parallel to the step edges. The wavelength of this 
wave pattern increases as the energy is lowered with re- 
spect to Ep. While a faint wave pattern discernable at 
-50 mV on terrace "3" (Fig. 0^) seems to originate from 
the subsurface impurity, the standing wave pattern at 
smaller energies is only perturbed by this impurity. The 
wave pattern does not originate from it. 

There is no obvious correspondence between wave pat- 
terns on adjacent terraces. For example, at -30 meV (Fig. 
^) the distance between the uppermost maxima is differ- 
ent on terrace "3" and on terrace "4". At -10 mV (Fig. 
0Jt) the wave patterns on terraces "3" and "4" are phase 
shifted in the upper half of the image. 

Fig. [S] compares the parallel wave pattern (-10 mV, 
Fig. 0J<) to a pattern at opposite polarity, i.e. in the 
unoccupied region (+90mV). At +90 mV the dominant 
wave pattern on terrace "3" runs perpendicular to the 
step edges, however, with contrast variations. Thus, per- 
pendicular and parallel wave patterns are superimposed. 
Note that even at +90 meV the narrower terraces, e.g. 
"5" , show the parallel wave pattern only. 

This superposition of parallel and perpendicular wave 
patterns becomes more pronounced on even broader ter- 
races (Fig.[5J). For the terrace "2" of 8.5 nm in width, the 
perpendicular wave pattern dominates at -38 mV (Fig. 
m>). At +15 mV and +30 mV the contrast variation 



along the perpendicular waves on the left half and on the 
right half of the terrace are independent of each other. 
This leads to an interference pattern at the middle of 
the terrace, which changes wavelength with energy (Fig. 
[5Jj to d). Note the differences of terrace "3" from the 
slightly narrower terrace "5" at + 15 mV and +30 mV. 
Terrace "5" shows lower conductivity at its center than 
at its sides, while terrace "2" shows higher conductivity 
at its center compared to its sides. Also terraces "3" and 
"4" show wave patterns of opposite contrast. 

The wavelength dependence on energy of the parallel 
wave pattern is investigated in more detail on the unoc- 
cupied side of the LDOS (Fig. 0). The standing wave 
pattern parallel to the surface steps on terrace "3" at 
positive polarity is contrast-reversed with respect to the 
standing wave pattern shown in Fig.0J The dl/dV map 
at -30 mV (Fig. |7p) shows that this is an effect of the 
opposite polarity of the tip-sample voltage bias and not 
of the terrace width. Due to a subsurface impurity (see 
below) the wave pattern in the lower part of this image 
is less regular. Thus, we concentrate on the upper part 
of the image. 

The average wavelength changes with bias voltage 
(Fig. 0ii) . For an ideal terrace with impenetrable edges, 
the isotropic free surface dispersion relation E(k) = 
E + 27^77 > with E o = -0.065 eV, is modified by "con- 
finement" due to the terrace structure, 

E n m)=Eo + - - + —JL, (1) 

where /c|| = 27r/A|| and n is the number of transverse 
modes. This equation predicts that for a fixed energy, 
narrower terraces with smaller width L will have a larger 
parallel wavelength, Ay. Examples of this effect can be 
seen in Figs.|l|to|7J although the oscillations in the nar- 
rower terraces are less regular for the narrower terrace. 
We attribute this to larger relative variations in width 
for fixed edge variations in the narrower terraces. 

Furthermore, a quadratic variation of A|| with E is ex- 
pected for a free electron dispersion (see Eq. Q) provided 
no additional transverse modes are occupied over the en- 
ergy range probed. Indeed, the data can be well fitted 
by a(E - Eo)- 1 / 2 . We find E = -(67 ± 5) meV, the 
same value as the bottom of the surface-state band de- 
termined from the surface-state onset on large terraces 
(see above) and from perpendicular wave patterns on 
large terraces^ Also the wavelengths determined from 
the negative dl/dV maps of Fig. 0]fit a quadratic dis- 
persion with E Q — —(68 ± 19) meV. Fig. [7Ji displays the 
dispersion relation deduced from this measurement that 
reproduces the dispersion relation determined from per- 
pendicular wave patterns^ Thus, parallel wave patterns 
originate also from the scattering of surface-state elec- 
trons. 
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III. THEORY 



A. Model and approach 



In this paper we are focused on understanding the role 
of disorder in the terrace edge profile on the observed 
modulations in the local density of states. Disorder in 
the edge profile breaks translational symmetry along the 
direction parallel to the terrace edges. This broken sym- 
metry prevents a simple separation of the energy (and the 
wavefunctions) of the electrons into a parallel and trans- 
verse component, as was done in Eq. In order to 
conveniently handle the broken translational symmetry 
another approach is needed. We chose to use a scatter- 
ing theory2& that has been quite successful in predicting 
the standing wave patterns in the LDOS observed with 
the STM around impurities and step edges on the sur- 
faces of the noble metals that support surface-states. 

The scattering theory given in Ref. |2f| most natu- 
rally describes the scattering of electrons in the surface- 
state from point-like impurities (adatoms) where an s- 
wave scattering approximation can be made; scattering 
from an extended impurity, such as a terrace edge, re- 
quires some modification of the theory, which otherwise 
cannot directly describe scattering from an extended ob- 
ject. Fortunately, since the standing wave patterns we 
are interested in (those parallel to the terrace edges) are 
observed most clearly in the central portion of the ter- 
race, we may make a simple approximation in the scat- 
tering theory for the terrace edge: We treat the terrace 
edge as a dense "row" of point (s-wave) scatterers. Far 
from the terrace edges (near the middle of the terrace, 
for example) the results will be indistinguishable from 
a terrace edge that is continuous. A "row" of atoms 
is expected to be a good approximation to the terrace 
edge provided the atoms that make up the wall are well 
within a wavelength (the wavelength of the surface-state 
electron) of each other and the distances probed are of 
the order of the separation of the atoms or larger. In the 
calculations described below our objective will be first to 
understand the qualitative features we can expect from 
scattering theory and second to compare as directly as 
possible with experiment. The latter requires determin- 
ing the phase shift of a single atom (one of many that 
make up the terrace edge) so that the correct phase shift 
results for the edge as a whole. Since the input to the 
scattering theory is the wall profile, the dispersion of the 
surface-state electrons, and the s-wave scattering phase 
shift 5, the only free parameter we have in the theory is 
5. We tune 8 for a fixed wall profile and surface-state 
dispersion to reproduce the experimental dl/dV vs. V 
spectra at the center of the terrace. Once S is optimized 
in this way, we computed the dl/dV image maps (at a 
fixed voltage) over all positions using this value and the 
exact edge profile of the terraces. 



B. Theoretical calculations 

1. Model calculation for qualitative effects 

To understand the qualitative effects we expect from 
a scattering theory approach to the LDOS modulations 
from disordered terraces, we first carried out calculations 
on a single terrace with perfectly straight walls and then 
added in increasing amounts of disorder. The terraces 
were modeled by placing 21 point (s-wave) scatterers at 
a distance of A_f/4 along the terrace edge. See Fig. [S] 

The dl/dV maps were computed using the surface- 
state dispersion E(k) = E + £ n k ff ■ (Confinement effects 
that do not cause an apparent shift in the surface-state 
band minimum are captured automatically by electron 
scattering from the terrace "edges".) An approximation 
to the "constant current" condition used in the exper- 
iment (while simultaneously aquiring topographic and 
dl/dV maps) is achieved by calculating 
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maps. 

In these initial model calculations, each scatterer is as- 
sumed to have a phase shift of tt/2. This value turns out 
to be almost "universal" for adatoms on the surfaces of 
noble metals2& so we take it as a starting point; later we 
will find a different value matches the experiments better 
when we attempt to make a direct, quantitative compar- 
ison between theory and experiment. The initial choice 
of 7r/2 is arbitrary and does not affect any of the qualita- 
tive conclusions reached from the calculations presented 
in Figs. IBland El While this value works well for single 
adatoms, the physics of scattering from a terrace edge 
is very different from the physics of scattering from an 
isolated atom, and this is why the best phase shift for 
the terrace edge atom differs from tt/2. 

Fig. |SJi shows the calculated [dI/dV] calc map for a 
trough of 5Xf in length and in width limited by 21 

point scatterers on each side. While the perpendicular- 
standing wave pattern outside the trough is clearly visi- 
ble, its narrowness inhibits a similar wave pattern within 
the trough. It is too small to support even half of a 
Fermi- wavelength. The introduction of an arbitrary de- 
fect in one of the terrace walls leads to standing waves 
parallel to the rows of atoms with the same wavelength 
Xp/2 as the perpendicular wave outside the trough (Fig. 
03}). Note the faint signs of oscillations in the "perfect" 
terrace as well. These originate from the finite length 
of the trough, where its ends act as defects. These first 
calculations demonstrate that a parallel wave pattern is 
indeed a result of terrace disorder. 

For a better understanding of particular features in the 
wave patterns observed in the experiments, we compare 
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[dI/dV] calc maps of varying terrace width and terrace 
disorder (Fig. EJ) ■ Terrace disorder has been modeled in 
the way described in the figure caption. Fig. a, d, and 
g support the idea that a wave pattern perpendicular to 
the steps can only develop for larger terraces. The inten- 
sity of the parallel wave pattern increases with increasing 
disorder (Figs. & to c) and shows the same qualitative 
features as in the experiments. For example, some of the 
maxima in Fig. |Sp are brighter than neighboring ones; 
some of them seem to be disconnected from neighboring 
maxima by deep minima, while between other maxima 
there are only little contrast variations. Fig. [5] c and i 
show that two bright protruding lines are seen instead of 
one for a larger terrace width. For the broadest disor- 
dered terrace (Fig. EJli, i) the two wave patterns are thus 
superimposed. Minima in the two perpendicular patterns 
are displaced with respect to each other by up to half a 
wave length (Fig. 0) as in the experiment. Many of the 
qualitative features of the experiments can thus be repro- 
duced by the scattering theory, showing they result from 
terrace edge disorder. 



2. Direct comparison with experiment 

The calculations presented so far were meant to illus- 
trate that the main physics is captured by the scatter- 
ing theory and to emphasize the quality of the quali- 
tative agreement between the scattering theory and ex- 
periment. Now, we go beyond these simple qualitative 
model calculations and use the actual surface-state dis- 
persion for Ag(lll) listed above and a terrace edge pro- 
file taken from experiment. The surface-state onset on a 
large terrace lies at E = —0.065 eV; the effective mass 
is m e ff = 0Am e . These values lead to Xf = 7.38 nm 
for the Fermi-wavelength on an infinite terrace. First, 
we calculate position dependent [dI/dV] calc spectra for a 
7.9 nm wide terrace with varying phase shift and compare 
them to the experimental spectra for the broadest terrace 
in Fig. The phase shift of the scatterers is chosen to 
best reproduce the Fabry-Perot type oscillations with the 
tip-sample bias, as shown in Fig. 1101 The best agreement 
is reached with a phase shift of 6 = (0.9 ± 0.1)7r. Those 
spectra resemble the experimental ones with a shift in the 
apparent surface-state onset and a second peak in the un- 
occupied region. Note that in order to keep our fitting 
parameter to only the number 5, we do not include any 
energy dependence in the phase shift. 

There is good agreement between the spectra calcu- 
lated for the right half of the terrace fFig. \Wh). if spec- 
trum number 14 showing no second maximum is taken 
as the central spectra, though it is measured 0.6 nm to 
the right of the geometrical center. This might be simply 
due to difficulties in determining the geometrical center 
of the terrace edges or may be due to the asymmetry 
of the reflection properties of ascending and descending 
step edges^ The reflectivities of steps on this surface 
have been measured to be r^esc — 0.72 — 2.91E/eV and 



fasc = 0.48 — 2.86E/eV— for the descending and the as- 
cending step, respectively. Spectra 15, 16, and 17 are 
measured at a distance of 0.6 nm, 1.2 nm, and 1.8 nm 
from spectra 14 are compared to calculations performed 
at '0.7 nm', '1.5 nm', and '2.0 nm' from the center of 
the terrace, respectively. In all cases, the peak at w 20 
meV lies at almost exactly the same energy and its rela- 
tive intensity is well reproduced. The first peak lies also 
at approximately the correct energy. However, a sharp 
rise in the theoretical spectra (an artifact of the calcula- 
tion procedure) makes a comparison difficult for the first 
peak. In both theory and experiment, the second peak is 
broader than the first peak, because the reflectivity and 
thus the lifetime of the resonance decreases with increas- 
ing energy. 

In Fig. 1101 the peaks are sharper and more intense in 
theory than in experiment. Several processes may ex- 
plain this difference in sharpness and intensity. On the 
experimental side, the spectra show thermal and modu- 
lation voltage broadening. We simulate the modulation 
voltage averaging and temperature broadening by aver- 
aging the theoretical spectra (Fig. HUb). Indeed, this 
broadening makes the additional peak in the theory at 
rj — 20 meV undetectable in the experiment and lowers 
the intensity of the peak at +20 meV. However, the theo- 
retical peaks remain somewhat sharper, presumably from 
details associated with how the terrace edges are mod- 
eled. One may think that the "porous" nature of the 
walls my lead to a shorter lifetime, and hence a broader 
thoery peak than what is observed experimentally. How- 
ever, this appears not to be the case. Changing the num- 
ber of atoms that make up the edge profile by 50% has a 
negligible effect on the widths of the peaks. 

Discrepancies may result from not correctly modeling 
the energy dependence of the reflectivity by replacing the 
steps by point scatters: The peaks in dl/dV can be asso- 
ciated with Fabry-Perot type resonance with an energy 
dependent reflectivity^ As a narrower resonance means 
a longer lifetime, and thus a greater reflectivity, we con- 
clude that theory has an effective reflectivity which is 
larger than experiment. The calculation has an energy 
dependent reflectivity, but because the walls are built 
up from point scatterers, it is difficult to determine the 
reflectivity exactly in order to compare it directly with 
experimental results for a "perfect" terrace. With effort, 
a very close match could be made by making the phase 
shift energy dependent and possibly adding imaginary 
parts as well, but extremely accurate numerical agree- 
ment is not our goal here. Rather, we wish to under- 
stand what is responsible for the essential features of the 
data. A more accurate calculation than ours would also 
include changes to the surface-state band dispersion it- 
self due to the terrace edges. Recall that our theory uses 
the surface-state dispersion of the free surface-states. 

The dl/dV spectra to the left and to the right of the 
center differ in both experiment and theory as a result 
of an asymetry in the disorder of the two terrace edges. 
However, for the experimental spectra to the left of the 



6 



center, the second peak shifts with position on the ter- 
race. This is not recovered in the calculation (Fig. HDb). 
We attribute this difference to the fact that the spectra 
are calculated for a terrace that is some 10% smaller than 
the experimental terrace and that has a different shape 
in detail, on the one hand, and to the larger difference 
in reflectivity properties between a row of atoms and a 
descending step than a row of atoms and an ascending 
step, on the other hand. The asymmetry in scattering 
properties of the two steps that can not be captured by 
replacing them by two rows of equivalent atoms leads 
thus to the 'off-center middle' (noted earlier in this sec- 
tion) and to larger discrepancies on the half of the ter- 
race that is closer to the descending step edge. Despite 
these deviations in detail, the scattering theory does sat- 
isfactorily reproduce the experimental spectra. Thus, we 
have calculated dl/dV maps for a particular terrace for 
a direct quantitative comparison of experiments to the- 
ory. We have chosen terrace "3" of Fig.0]due to the clear 
wave pattern observed. This terrace is also broad enough 
that changes to the surface-state on narrow terraces^ are 
negligible. 

The terrace walls are approximated by 54 atoms placed 
as close to the actual boundary as possible (Fig. Illh . 
white dots). Fig. EU compares the [dI/dV] mcas maps, 
computed from Eq. Q), of Fig. 0] and k at -20 mV, at 
-10 mV, and of Fig. at +90 mV to the calculation. 
The qualitative agreement is excellent: A similar wave 
pattern is observed in the calculated spectra as in the 
experimental one, both in the occupied and the unoccu- 
pied region. Also the decrease in wavelength from -20 
meV to -10 meV is reproduced in the calculation. The 
comparison of the [d!/dV] mcas line scans through the cen- 
ter of the terraces (Fig. II lb ) shows this good wavelength 
agreement for both -20 meV and -10 meV. The fourth 
maximum at -20 mV in Fig. II lb is simply not resolved in 
the experimental data. Only the relative intensities are 
partly recovered. 

We have also simulated the broadest terrace "3" in Fig. 
for a more detailed comparison of theory with exper- 
iment for wavelength changes in the unoccupied region. 
Fig. 1121 shows the comparison of simulation and exper- 
iment at 33, 26, and 17 mV. Unfortunately, the direct 
comparison is somewhat obstructed by an impurity, vis- 
ible by enhancing the contrast in the topographic image 
(Fig. \V2h ) leading to a very particular wave pattern at 
-30 meV (Fig. |7p). At first sight, the comparison seems 
to be less favourable. The simulation shows two bright 
protruding lines close to the step edges separated by a 
dark stripe with little contrast variations. However, this 
reproduces the variation of intensity across the terrace in 
the experiments (from bright to weaker to bright). De- 
spite the impurity the lines cans through the middle of 
the terrace agree well in wavelength in the upper part of 



the image fFig.11211'). Only a phase shift in the lower part 
is not recovered. Thus, again the major difference lies in 
intensity. This last example underlines the complexity of 
the evolved wave pattern and that it may depend quite 
sensitively on various details of the geometry. 



IV. CONCLUSIONS 

We have probed the LDOS of disordered terraces of 
varying width and observed oscillations parallel to the 
terraces with a wavelength in accordance with a free elec- 
tron dispersion, and we have shown that these parallel os- 
cillations arise naturally in scattering theory calculations 
for disordered terraces. The theory also captures impor- 
tant contrast variations with tip-sample voltage polarity. 
Thus, the wave patterns observed experimentally are a 
result of step disorder. 

The complete pattern is a combination of Fabry-Perot 
(surface-state) and disorder scattering effects along with 
surface reconstruction effects. The Fabry-Perot effects 
are intrinsic in scattering theory, but the changes to 
the underlying surface-state band structure are not (in 
the theortical model used here, which assumes a clean 
surface-state dispersion). As the scattering theory re- 
covers all the qualitative features of the experiment 
and yields semi-quantitative agreement, the surface-state 
scattering effects dominate (at least for the wider ter- 
races). Remaining differences are attributed to surface 
reconstruction effects. The semi-quantitative agreement 
in the simulation also shows that the cross-talk between 
the large terraces and their neighbors is negligible. 

The terraces that the dl/dV maps have been calcu- 
lated for are wide enough to allow us to use the surface- 
state band structure valid for infinitely broad terraces. In 
principle narrower terraces could be calculated within the 
same theory by using the appropriate dispersion relation 
for narrower terraces, which will in general be anisotropic 
and depend on terrace widths 

As a future project, we suggest to locally investigate 
the reaction dependence on the vicinal surface and we 
expect local variation on the nanometer scale due to the 
reported variation of the LDOS. 
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FIG. 1: (a) Real space STM image taken at 7 K showing spa- 
tial oscillations of electrons in the Ag(lll) surface-state on 
a large terrace 51 nm in width. Circularly symmetric oscil- 
lations in the LDOS are seen around point impurities. The 
tip-sample voltage bias is Ubi as =32 mV with a tunneling cur- 
rent hun =0.43 nA. (b) The surface state onset is indicated in 
the dl/dV spectrum on the large terrace shown in (a). Here 
Von = {Utop + Ubot)/2 is taken as the value of the onset, (c) 
The dl/dV spectrum on the edge of a narrow terrace of 6.1 
nm width (indicated by a white dot on the inset STM image) 
shows two peaks in intensity. These peaks can be related to 
the onset of the lowest transverse mode of the terrace (first 
peak fa —50 mV) and the first excited tranverse mode (second 
peak « 50 mV). Here [4i as =-137 mV and J tun =0.34 nA. 




FIG. 2: (Color online.) Position dependent STS at differ- 
ent points across the width of different terraces showing the 
dependence of the surface-state onset and of the second max- 
imum in dl/dV spectra on average terrace width and on po- 
sition on the terrace, (a) The background grayscale image 
is a topographic image (?7(,j aa =-47 mV, / ttl „=1.3 nA) show- 
ing approximately 5 terraces with irregular, disordered edges. 
The average widths range from 2.7-7.9 nm. Superimposed 
on each terrace of the STM image is a set of inset STS data 
taken at various positions (indicated in color) across the ter- 
race. The positions are numbered from left (1) to right (30). 
(b) Position of bottom (down triangle), top (up triangle), and 
middle (diamond) of first (filled symbols) and second (open 
symbols) peak in dl/dV spectra of (a). The second peak is 
only observed on the widest terrace (of width 7.9 nm) and is 
related to exciting the second transverse mode. Dark vertical 
lines indicate positions of the step (terrace) edges, (c) Rela- 
tive intensity [i.e. (dl / dVt op — dl/ dVbot) /dl / dVbot] of the first 
peak (filled symbols) and the second peak (open symbols) for 
spectra of (a). 
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FIG. 3: Scaling of the relative position of the first and sec- 
ond maximum in dl/dV curves on small terraces vs. inverse 
terrace width squared. For terraces wider than L ~ 4.3 nm, 
\U° n — [/J™ | L~ 2 as shown by the linear fit, implying an in- 
terpretation consistent with the presence of transverse modes 
of the surface state. However, the effective mass determined, 
m e77 = (0-59 ± 0.17) m e , deviates from the effective mass of 
electrons in the surface-state band where m e // = 0.4m e . On 
the other hand, terraces narrower than L ~ 3.4 nm do not 
scale linearly with L~ 2 , but instead show behavior consistent 
with step modulation. 
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FIG. 4: dl/dV maps of Ag(lll) terraces in the occupied re- 
gion of the spectrum [negative tip-sample voltage bias, in- 
dicated in the upper left-hand corner of (b)-(k)]. The AC 
voltage modulation is U mo d = 4 mV and the frequency is 
Vmod = 738.1 Hz. (a) Topographic image (C/f,i as =-100 mV, 
Itun=l3 nA). The two central terraces are labeled "3" and 
"4". The arrow points to an impurity at step a edge be- 
tween the two. (b)-(k) Show dl/dV maps for Ubias = -100 
mV to -10 mV as indicated. Note the difference in wave 
patterns on the different terraces. In particular, the ter- 
races "3" and "4" begin to develop oscillations in the LDOS 
parallel to the step edges over the range of tip-sample bias 
—40 meV < Ubias i$ — 10 meV. The oscillations on terrace 
"3" are investigated in more detail in Fig. 1111 



FIG. 5: Comparison of dl/dV maps for different tip-sample 
polarity with AC modulation parameters U mo d = 4 mV, 
Vmod = 738.1 Hz and Itun = 1.4 nA. (a) Topographic image 
of approximately 6 irregular terraces. The terraces are the 
same as those shown in Fig. [I] (upper half), (b) dl/dV map 
at U hias = -10 mV. (c) dl/dV map at U bias = 90 mV. Note 
the parallel and perpendicular wave pattern on terrace "3" 
for negative and positive polarity, respectively. This change 
in the wave pattern can be understood with the scattering 
theory described in the text. See Fig. ^3 




FIG. 6: Evolution of terrace wave patterns with Ubias- The 
terraces shown are different from the set shown in Figs. |3] 
and|^| The STM parameters are U mo d = 4 mV, v mo d = 327.9 
Hz and itun=l nA. (a) Topographic image taken at Ubias = 
15 mV, it„„=l nA. Terraces "2", "3", "4", and "5" develop 
distinct wave patterns with varying Ubias, indicating a clear 
dependence of wave pattern on both terrace edge disorder and 
terrace width, (b)-(d) dl/dV maps at -38 mV, +15 mV, and 
+30 mV. Note the interesting superposition of parallel and 
perpendicular wave patterns on terrace "2". 
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FIG. 7: Wavelength dependence of parallel wave patterns in 
unoccupied region of the energy spectrum (positive tip-sample 
bias). Compare with Fig. [I] which shows different terraces 
at negative Utias- STM parameters are Umod = 8 mV and 
v m od = 381.8 Hz. (a) Topographic image (-30 mV, 0.39 nA). 
(b-1) dl/dV maps for U bias = 17 mV to 100 mV and -30 mV 
as indicated. Black stripes (in h-1) are an artifact of the STM 
electronics and should be ignored. For better visibility (h)-(l) 
have four times the intensity contrast of (c)-(g). Note that 
there is a "contrast" reversal of the wave pattern of terrace 
"3" compared to those image shown in Fig. 0] This feature 
is also captured by the scattering theory. (See Fig. 1121 1 (m) 
Wavelength of standing waves on terrace "3" in (b) to (1) 
with square root fit. (See Sec. IIB in paragraph just below 
Eq. Q.) (n) Dispersion relation as deduced from the fit in 
(m). 




FIG. 8: Model scattering theory calculations illustrating cre- 
ation of standing wave patterns in the LDOS parallel to the 
terrace edges when disorder is present. Here "?" = "A" due 
to a graphics error. The terrace edges are modeled by 
rows of atoms. As long as the atoms are much closer together 
than the wavelength of the electrons in the surface state this 
is a good approximation. Shown are two [dl /dV] calc maps 
calculated with scattering theory for a terrace of length 5Xf 
and width Af/4 with (a) no defect and (b) one defect. The 
introduction of a single point defect creates a strong varia- 
tion in the LDOS parallel to the edges of the terrace. In both 
cases there are oscillations clearly visible perpendicular to the 
terraces edges outside the terrace. 
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FIG. 9: Theoretical dl/dV maps for different width and dis- 
order. Disorder is introduced manually by shifting the hor- 
izontal positions of the atoms that model the walls of the 
terrace. The characteristics of the disorder (such as whether 
or not the displacements are Guassian distributed, what the 
RMS value of the fluctuations is, etc.) are unimportant for 
the qualitative study here. The maximum displacement of 
the "random" positions is: (a,d,g) 0, i.e., no disorder, (b,e,h) 
±8% variation in width, and (c,f,i) ±20% variation in width. 
The average widths of the terraces are: (a,b,c) 0.8Af, (d,e,f) 
Af, and (g,h,i) 1.2Af- The features seen in (c) and (i) are 
qualitatively similar to those seen on terrace "3" in Fig. 0] at 
negative and positive voltage biases, respectively. See Fig. 1111 
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FIG. 10: (Color online.) Comparison of dl/dV spectra taken 
at different positions in the second terrace in Fig. [5^ of width 
7.9 nm. Solid lines: calculation with phase shift 0.97T in the 
surface-state scattering theory described in the text. Dotted 
lines: Experimental results. (a,c) Scattering calculation with- 
out thermal broadening, (b) same as (a), but theoretical curve 
is averaged over an energy window to simulate broadening due 
to thermal effects and modulation voltage. The quality of the 
agreement between experiment and theory in (b) is good con- 
sidering the range of energy and the relative simplicity of the 
theory. The main disagreement occurs around U w — 50 meV 
where the theory predicts a narrower and higher peak than is 
seen experimentally. This indicates the approximations used 
to model the terrace edge result in a more highly reflective 
edge at these energies than the actual edge in experiment. 
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FIG. 11: Comparison of theoretical to experimental dl/dV 
maps for terrace "3" of Fig. 2] f° r bias voltages as indicated 
above the images, (a) Theoretical dl/dV maps computed us- 
ing Eq. (b) Experimental dl/dV maps. Note that the 
field of view in y-direction is larger in the theoretical maps 
than in the experimental maps as indicated in (a) by the faint 
horizontal lines near the top and bottom of the figures. The 
region in between these horizontal lines should be compared 
with the experimental figures in (b). (c) We show line scans 
through the middle of terrace "3". Gray: theory. Black: ex- 
periment. The y-axis is different in the two spectra for better 
comparison. 
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FIG. 12: Comparison of theoretical to experimental dl/dV 
map for unoccupied states (positive tip-sample bias), (a) To- 
pography with enhanced contrast. Black dot on terrace in- 
dicates a point impurity, (b) Calculated dl/dV maps. As 
in Fig. 1111 the region between the two faint horizontal lines 
should be compared with the experimental data in (c). (c) 
Experimental dl/dV maps of Fig. Q (d) Line scans through 
middle of terrace. Gray: theory, Black: experiment. 



